NAVAL  POSTGRADUATE  SCHOOL 
MONTEREY,  CALIFORNIA 


THESIS 

SIMPLIFIED  MICROMECHANICAL  MODELS  FOR 
ANALYSIS  OF  INTERFACE  DEBONDING  IN  A 
FIBROUS  COMPOSITE 

by 

Guo,  Jia-Yuam 
December,  1995 

Thesis  Advisor:  Young  W.  Kwon 

Approved  for  public  release;  distribution  is  unlimited. 


1 996041 5  079 


REPORT  DOCUMENTATION  PAGE 


Form  Approved  0MB  No.  0704-0188 


Public  reporting  burden  for  this  collection  of  information  Is  estimated  to  average  1  hour  per  response,  including  the  time  for  reviewing  instruction,  searching 
existing  data  sources,  gathering  and  maintaining  the  data  needed,  and  completing  and  reviewing  the  collection  of  information.  Send  comments  regarding 
this  burden  estimate  or  any  other  aspect  of  this  collection  of  information,  including  suggestions  for  reducing  this  burden,  to  Washington  Headquarters 
Services.  Directorate  for  Information  Operations  and  Reports,  1215  Jefferson  Davis  Highway,  Suite  1204,  Arlington,  VA  22202-4302,  and  to  the  Office  of 
Management  and  Budget,  Paperwork  Reduction  Project  (0704-0188)  Washington  DC  20503.  _  _ 


1.  AGENCY  USE  ONLY  rLeave  2.  REPORT  DATE  3.  REPORT  TYPE  AND  DATES  COVERED 

blank)  December  1995  Master’s  Thesis 


4.  TITLE  AND  SUBTITLE  SIMPLIFIED  MICROMECHANICAL  MODELS  FOR  5.  FUNDING  NUMBERS 
ANALYSIS  OF  INTERFACE  DEBONDING  IN  A  FIBROUS  COMPOSITE 
6.  AUTHOR(S)  Guo,  Jia-Yuarn 


7.  PERFORMING  ORGANIZATION  NAME(S)  AND  ADDRESS(ES) 
Naval  Postgraduate  School 
Monterey  CA  93943-5000 


9.  SPONSORING/MONITORING  AGENCY  NAME(S)  AND  ADDRESS(ES) 


8.  PERFORMING 
ORGANIZATION 
REPORT  NUMBER 


10.  SPONSORING/MONITORING 
AGENCY  REPORT  NUMBER 


11.  SUPPLEMENTARY  NOTES  The  views  expressed  in  this  thesis  are  those  of  the  author  and  do  not  reflect 
the  official  policy  or  position  of  the  Department  of  Defense  or  the  U.S.  Government. _ 


12a.  DISTRIBUTION/AVAILABILITY  STATEMENT  12b.  DISTRIBUTION  CODE 

Approved  for  public  release;  distribution  is  unlimited. 


13.  ABSTRACT  (maximum  200  words) 

The  objective  of  this  study  is  to  develop  simplified  micromechanical  models  to  analyze  the 
interface  debonding  between  fiber  and  the  matrix  materials.  Both  analytical  and  simplified  finite 
element  models  are  used  to  predict  the  effective  transverse  elastic  modules  of  fibrous  composites 
with  a  partial  interface  crack  based  on  the  material  properties  of  their  constituents.  The  simplified 
finite  element  model  uses  springs  in  the  connecting  nodes  between  the  fiber  and  matrix.  A  detailed 
finite  element  analysis,  which  is  programmed  using  the  MATLAB  engineering  software  is 
performed  to  check  the  accuracy  of  the  simplified  models.  The  simplified  models  yield  accurate 
effective  transverse  elastic  modulii  of  various  composites  with  partial  interface  cracks  when 
compared  to  the  results  obtained  from  detailed  finite  element  analyses. 


14.  SUBJECT  TERMS  Micromechanical  Models,  Interface  Debonding.  Fibrous  Composite 


17.  SECURITY  CLASSIFI¬ 
CATION  OF  REPORT 
Unclassified 


18.  SECURITY  CLASSIFI¬ 
CATION  OF  THIS 
PAGE 

Unclassified 


19.  SECURITY  CLASSIFI¬ 
CATION  OF 
ABSTRACT 
Unclassified 


15.  NUMBER  OF 
PAGES  116 


16.  PRICE  CODE 


20.  LIMITATION  OF 
ABSTRACT 
UL 


NSN  7540-1-280-5500  Standard  Form  298  (Rev.  2-89) 

Prescribed  by  ANSI  Std.  239-18  298-102 


1 


11 


Approved  for  public  release;  distribution  is  unlimited. 

SIMPLIFIED  MICROMECHANICAL  MODELS  FOR  ANALYSIS  OF 
INTERFACE  DEBONDING  IN  A  FIBROUS  COMPOSITE 


Guo,  Jia-Yuam 

Major,  Republic  of  China  Army 
B.S.,  Chung  Cheng  Institute  of  Technology  - 1986 

Submitted  in  partial  fulfillment 
of  the  requirements  for  the  degree  of 

MASTER  OF  SCIENCE  IN  ENGINEERING  SCIENCE 

from  the 

NAVAL  POSTGRADUATE  SCHOOL 
December  1995 


Author: 


Approved  by: 


C/  <you 


bung  W.  Kwon,  Thesis  Advisor 


(  Matthew  D.  Kelleher,  Chairman 
Mechanical  Engineering  Department 


111 


iv 


ABSTRACT 


The  objective  of  this  study  is  to  develop  simplified  micromechanical  models  to 
analyze  the  interface  debonding  between  fiber  and  the  matrix  materials.  Both  analytical 
and  simplified  finite  element  models  are  used  to  predict  the  effective  transverse  elastic 
modules  of  fibrous  composites  with  a  partial  interface  crack  based  on  the  material 
properties  of  their  constituents.  The  simplified  finite  element  model  uses  springs  in  the 
connecting  nodes  between  the  fiber  and  matrix.  A  detailed  finite  element  analysis,  which 
is  programmed  using  the  MATLAB  engineering  software  is  performed  to  check  the 
accuracy  of  the  simplified  models.  The  simplified  models  yield  accurate  effective 
transverse  elastic  modulii  of  various  composites  with  partial  interface  cracks  when 
compared  to  the  results  obtained  from  detailed  finite  element  analyses. 
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1.  INTRODUCTION 


Use  of  fibrous  composite  material  is  expanding  rapidly,  especially  in  structural 
application  because  properties  of  fibrous  composite  materials  are  beneficial  such  as  high 
ratios  of  stiffhess-to-weight  and  strength-to-weight. 

Prediction  of  effective  modulii  of  fibrous  composites  from  material  properties  of 
their  constituents  has  been  a  main  topic,  especially  when  there  are  micro-cracks.  There 
are  many  different  micromechanics  models  which  have  been  proposed  to  predict  effective 
moduli.  However,  most  of  them  did  deal  with  neither  damages  nor  3-D  models. 
Broutman  and  Sahn  observed  the  transverse  matrix  cracking  in  cross-ply  composites 
[Ref.  1].  After  the  observation,  the  phenomenon  has  been  studied  anal5ftically  and 
experimentally  by  several  investigators.  Highsmith  and  Reifsnider  [Ref  2]  conducted  an 
experiment  to  measure  the  longitudinal  stiffness  reductions  of  various  laminates.  Talreja 
[Ref.  3,4]  used  continuum  damage  mechanics  to  represent  damage  as  internal  variables  in 
a  phenomenological  approach,  and  Kwon  [Ref  5]  used  an  analytical  micromechanics 
model  to  predict  the  effective  modulii  of  composites.  Kwon  and  Bemer  [Ref.  6]  also 
used  a  micromechanics  model  to  investigate  matrix  cracking  in  laminated  fibrous 
composite,  and  predict  matrix  cracking  and  its  effect  on  stiffhess  reduction. 

An  analysis  of  a  fibrous  composite  structure  requires  information  of  material 
properties  of  the  composite.  To  determine  the  material  properties  of  a  composite,  it  is 
necessary  to  conduct  experiments  that  will  allow  measurement  of  those  properties.  Even 
though  some  of  the  material  properties  can  be  determined  with  physical  tests,  it  is 
difficult  to  obtain  some  other  properties  of  the  composite. 

This  study  will  develop  simplified  micromechanical  models  to  analyze  the 
interface  debonding  between  the  fiber  and  the  matrix  materials  using  an  analytical  model 
and  a  simple  two-dimensional  finite  element  model  with  which  the  effective  transverse 
elastic  modulii  of  fibrous  composites  with  partial  interface  cracks  can  be  predicted.  Since 
modeling  an  interface  crack  in  the  simplified  finite  element  model  is  limited  by  a 
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minimum  number  of  node  points,  the  model  uses  springs  in  the  connecting  nodes 
between  the  fiber  and  matrix.  In  order  to  check  the  accuracy  of  the  two  simplified 
micromechanics  models,  a  detailed  finite  element  analysis  is  performed. 

Using  properties  of  the  fiber  and  matrix  materials  as  well  as  their  fiber  volume 
firactions,  it  can  be  determined  the  relationship  between  the  elastic  modulus  and  the  crack 
length.  Figure  1  shows  a  micromechanical  modeling  process.  In  later  chapters,  both 
analytical  and  finite  element  micromechanical  models  will  be  presented,  and  followed  by 
the  results  and  discussion.  The  last  chapter  has  conclusions  and  recommendations. 
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Figure  1.  Modeling  Process  of  a  Micromechanical  Model  From  a  Fibrous 

Composite. 


II.  MICROMECHANICS  MODELS 


A.  ANALYTICAL  MICROMECHANICS  MODEL 

This  micromechanics  model,  which  considers  a  unit  cell  consisting  of  the  fiber 
and  surrounding  matrix  material,  was  developed  by  Kwon  [Ref.  5],  and  is  shown  in 
Figure  2.  The  cross-section  of  the  fiber  is  assumed  to  be  square,  for  simplicity,  and  the 
unit  cell  is  divided  into  four  subcells  a,b,c  and  d.  Subcell  ‘n’  represents  the  fiber  section 
and  other  subcells  are  the  matrix.  The  size  of  the  unit  cell  is  unity,  and  the  size  of  the 
subcell  varies  with  the  fiber  volume  fraction.  The  stress  continuity  at  the  subcell 
interfaces  is  expressed  as: 

a  b  —C  _«/  h  _  _i/ 

^22  ^22  ’  ^22  “  ^22  ’  ^33  ^33  ’  ^33  ”  ^33 

^12  ”  ^12 ’^12  ""  ^12  ’  ^13  ”■  ^13’ ^13  ^13 

a  _  b  _  a  _  b 
^23  ““  ^23  ^  ^23  ^23  ’ 

and  the  strain  compatibility  is  assumed  as: 

£u  = 

S22  +  £22  ~  ^22  ^22  ’  ^33  ^33  ~  ^33  ^33  (2*2) 

f "  +  Sn  =  £u  +  £n  >  ^13  +  ■ 

The  constitutive  relation  within  each  subcell  can  be  expressed  as: 


_cr  __  T^a  a 
^  if  —  ^ijkl^kl  5 


where  i,j,k,l  =  1,2,3  and  a  =  a,b,c,d. 

From  Figure  2  (before  crack),  Kwon  [Ref  5]expressed  the  composite  stresses  and 
strains  as  a  function  of  the  subcell  stresses,  subcell  strains  and  the  fiber  voltune  fraction 


as  follows: 
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where  and  f,^are  composite  stresses  and  strains,  cr"//  and  f“;/  are  stresses  and  strains 

of  each  subcell,  and  Vf  is  the  fiber  volume  fraction  that  determines  the  size  of  the  subcell. 
The  composite  stresses  and  strains  are  determined  from  the  volume  average  of  the 
stresses  and  strains  of  the  subcells. 
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a  -  Fiber  Subcell 

b,  c  and  d  -  Matrix  Subcells 


Figure  2.  Subcells  of  a  Unit  Cell. 


When  there  is  an  interface  crack  between  the  fiber  and  the  matrix  as  shown  in 
Figure  3,  some  volume  portions  of  subcells  a  and  c  cannot  sustain  a  load  applied  in  the  3- 
direction.  As  a  result,  equation  (2-4)  is  rewritten  as 


6 


(2-6) 
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a  -  Fiber  Subcell 

b,  c  and  d  -  Matrix  Subcells  t  Crack  Length 


Figure  3.  Crack  Situation  in  Unit  Cell. 


In  order  to  find  the  relation  between  composite  stresses  and  composite  strains, 
first  of  all  equations  (2-1)  through  (2-3)  and  (2-5)  are  solved  together.  Then,  the 
following  matrix  equation  can  be  obtained. 

=  =  1,2,3  ,  (2-7) 

where  [A]  is  a  constant  coefficient  matrix  and  is  invertible. 

Pre-multiplying  [A] on  both  sides  of  equation  2-7,  we  get: 
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Finally,  combining  equation  2-4  and  equation  2-8  along  with  equation  2-3  yields  the 
relationship  between  composite  stresses  and  strains: 

{5„}=[C]{s„}  .  (2-9) 

where  [C]  is  the  smeared  material  property  matrix  of  the  eomposite.  In  case  of  the 
interface  crack,  equation  (2-6)  is  used  instead  of  equation  (2-4)  with  the  same  procedures. 

In  Kwon's  model,  he  assumed  the  fiber  and  matrix  were  orthotropic  materials. 
This  means  the  normal  stress  (strain)  components  and  the  shear  stress  (strain)  components 
are  uncoupled.  Appendix  A  shows  the  equations  relating  the  subcell  strains  to  the 
composite  strains  for  the  normal  strains. 

B.  FINITE  ELEMENT  MICROMECHANICS  MODEL 

The  finite  element  method  (FEM)  is  a  numerical  procedure  for  solving  differential 
equations  of  physics  and  engineering.  A  continuous  quantity  such  as  displacements, 
temperature,  or  pressure  can  be  approximated  by  a  piece-wise  continuous  function 
spanned  over  a  finite  number  of  subdomains  into  which  the  problem  domain  is  divided. 
For  each  subdomain,  nodal  variables  are  defined  as  unknowns  which  satisfy  the 
governing  equation  in  an  approximate  manner.  The  nodal  unknowns  are  assembled  into  a 
system  of  matrix  equation  for  the  whole  domain.  The  matrix  equation  is  solved  for  the 
unknowns  with  prescribed  boundary  conditions. 

1.  Finite  Element  Method 

In  this  study,  one  of  the  finite  element  methods  called  the  weighted  residual 
method  (WRM)  has  been  selected  for  the  micromechanical  model.  The  concept  of  this 
method  is  to  minimize  the  residual  multiplied  by  a  weighting  function  (W).  The  WRM 
has  many  forms,  and  one  that  is  very  useful  for  solving  boundary  value  problems  is 
Galerkin's  method.  This  is  the  method  used  in  this  study.  The  proeedure  for  this  method 
is  given  below: 

The  equilibrium  equations  for  two  dimensional  elasticity  are: 
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^<^x 

=  0  , 

dy 

dx 

yx 

+ 

=  0  , 

dx. 

4^ 

(2-10) 


and  the  shear  stresses  satisfy: 


^  xy  ^  yx  »  (2“11) 

The  two-dimensional  strain-displacement  relationships  i.e.  (kinemetic  equations)  are: 
du  dv  du  dv 

=  r,.  =^  +  —  >  (2-12) 


du 


_du  3^ 


where  u  and  //  represent  the  displacements  in  x  and  y  direction. 

The  stress-strain  relationships  i.e.  (constitutive  equation)  for  the  plane  stress  condition 
are: 


- - ^  V  1  0 

1-//  \-u 

0  0  y 

2  Vxy 


(2-13a) 


Where  E  and  //  are  elastic  modulus  and  Poisson’s  ratio. 

The  traction  boundary  conditions  are: 

<l>x  =  (^x^x  +  ^xy^y 

<l>y  =  '^xynx^<^y^y  ' 

The  basic  formulation  of  the  weighted  residual  method  is: 

l=\RW^d^  =  0  (/  =  U)  , 

n 

where  R  is  residual  of  the  governing  equation,  and  Wi  is  a  weight  function. 
By  applying  equation  2-10  to  the  weighted  residual  equation  we  get: 


(2-13b) 


(2-14) 


(2-15) 
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Using  integrations  by  parts  and  combining  the  two  I  functions  yields: 


- -  +  ^xv 

r  3c  ^  ^ 

J  3(V,  3fV, 

3c  ^  ^ 


m  1 


(2-16) 


(2-17) 


0  ^ 


n  0  ^  ^  \r^ 

dy  dx. 


rn\m 


,  (2-18) 


where  <|)x  and  are  boxmdary  tractions  as  stated  in  equation  2-14.  Here,  /2  denotes  the 
domain  and  /’indicates  the  boundary. 

Plugging  the  strain-displacement  equations  into  the  constitutive  equations  and 
substituting  the  resultant  equation  into  2-18  yields: 


j  3c 
n  0 


3fV^  31V, 
^  3c 


1^1  “  It 

Z  1. 

dy  3c 


(2-19) 


In  this  model,  linear  triangular  elements  are  used  for  the  formulation.  That  is,  each 
triangular  element  has  three  nodes  and  each  node  has  two  displacements  u,  v  in  the  x  and 
y  direction  so  that  each  element  has  six  degrees  of  freedom.  The  shape  function  that  is 


10 


used  for  the  triangular  element  is  Hi  (x,y).  The  displacements  are  interpolated  using  the 
shape  functions. 


or  in  matrix  form: 


u^±H,u, 

/=I 


(2-20) 


0  0  H,  0 

0  H  0  K  0  H 


(2-21) 


The  displacement  vector  is  2  x  1,  cmd  can  be  expressed  as  the  shape  function 
matrix  (2x6)  times  the  nodal  displacement  vector  (6x1);  as  shown  above.  We  define  a 
matrix  [B]  that  is  obtained  by  taking  first  partial  derivatives  of  the  shape  functions  with 
respect  to  the  coordinates,  that  is: 


3c 


0 

0 

3H, 

0 

3c 

3c 

3H, 

0 

3H^ 

0 

cH, 

3H^ 

3H, 

3H, 

3c 

3c 

3k 

(2-22) 


Because  the  weighting  function  (Wi)  of  the  Galerkin  method  is  the  same  as  the  shape 
function: 


(2-23) 


the  weighting  fimction  matrix  in  equation  2-19  can  be  written  as: 
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or 

[A:,]{rf)=[iinD][fi]A(4  =  {F4  ,  (2-27) 

where  A  is  the  area  of  the  triangular  element,  t  is  the  thickness  of  the  element.  For  unit 
thickness,  t  is  equal  to  1 . 

The  right  hand  side  of  the  equation  2-27  is  called  the  element  force  vector  {Fg} 
which  is  dependent  on  the  surface  traction  boundary  conditions  and  represents  the  total 
force  at  the  nodal  points  in  the  x  and  y  directions.  Combining  the  element  stiffiiess 
matrices  into  the  system  global  matrix  equation  gives: 
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[KK]  [dd]  ={FF} 


(2-28) 


The  FEM  used  in  this  study  is  to  calculate  the  elastic  modulus  of  the  composite 
material,  whether  there  is  a  crack  or  not.  From  equation  (2-28)  we  can  solve  for  the 
nodals  displacement  {dd}  by  applying  the  displacement  boundary  condition.  After  we  get 
tfie  nodals  displacements,  stresses  are  computed  for  each  element.  The  elastic  modulus 
can  be  calculated  from: 

.  (2-29) 

where  is  the  composite  average  stress  and  e  ^  is  the  composite  average  strain.  In  this 
study,  we  applied  a  uniform  displacement  to  the  micromechanics  model  as  shown  in 
Figure  4  as  the  displacement  boundary  condition. 
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a  -  Fiber  Subcell 

b,  c  and  d  -  Matrix  Subcells 


Figure  4.  Unit  Displacement  Model. 

Two  different  kinds  of  finite  element  analyses  are  performed.  The  first  one  is  a 
simplified  model  in  which  each  subcell  in  Figure  4  is  divided  into  two  triangular 
elements.  The  other  one  is  a  detailed  model  for  which  each  subcell  is  divided  into  many 
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niimber  of  finite  elements.  The  latter  model  was  used  to  check  the  accuracy  of  the 
simplified  analytical  and  finite  element  models. 


2.  Interface  Crack 

When  a  partial  interface  crack  between  the  fiber  and  the  matrix  exists,  the  crack 
caimot  be  modeled  explicitly  for  the  simplified  finite  element  model  because  there  is  no 
node  to  represent  the  crack  tip.  As  a  result,  it  is  assumed  that  there  are  elastic  springs  at 
the  two  comer  nodes,  which  connect  the  fiber  and  the  matrix.  The  double  node  technique 
is  used  for  this  purpose.  Since  the  nodes  have  two  degrees  of  freedom  (x  and  y 
directions),  the  spring  constants  at  every  node  are  kx  and  ky.  The  stiffiiess  matrix  for  the 
two  dimensional  spring  element  is; 


\  0  -k,  0 

0  A,  0  ky 

-k.  0  A:,  0 


Q  -k  Q  ky 


(2-30) 


By  adding  the  stiffness  matrix  into  the  global  system  matrix,  springs  connect  the 
fiber  to  the  matrix.  The  spring  constant  determines  the  status  of  a  partial  interface  crack. 
For  example,  if  the  spring  constant  is  very  large,  there  is  no  crack  and  if  the  spring 
constant  is  zero,  there  is  a  complete  crack  at  the  interface.  An  intermediate  value  of  spring 
constant  represents  a  partial  crack.  The  proper  spring  constants  need  to  be  determined 
later  for  a  given  size  of  interface  crack. 

C.  MATLAB  COMPUTER  PROGRAM  DEVELOPMENT 


The  MATLAB  engineering  software  is  chosen  to  program  the  finite  element 
procedure,  and  to  calculate  the  composite  material  property  of  a  composite  with  perfect 
bonding  or  partial  cracks  between  the  fiber  and  matrix.  The  following  procedures  are  the 
outline  of  the  MATLAB  program  development: 

1 .  Define  the  numbers  of  elements  and  nodes  and  input  x,  y  coordinates  of  the  nodes 
associated  with  each  element  using  'GRID'. 
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2.  Input  the  material  properties  {E  &  v)  of  the  fiber  and  matrix,  and  calculate  the 
kinematic  matrix  [B],  the  material  property  matrix  [D]  and  the  element  stiffness 
matrix  [K^]  using  ‘TESTT’. 

3.  Use  ‘CHECK’  to  assemble  each  element  stiffness  matrix  [K^]  into  the  system 
stiffness  matrix  [KK]. 

4.  Use  ‘MODTT’  to  apply  the  geometric  boundary  conditions  and  to  modify  the 
system  global  stiffness  matrix  [KK]. 

5.  Compute  the  displacement  {dd}  and  stresses  using  “STRES.” 

6.  From  the  main  program  “MAIN,”  we  obtain  the  composite  material  property,  and 
compare  the  results  by  plotting. 

7.  For  the  interface  crack  model,  we  add  the  spring  stiffness  matrix  to  the  system 
stiffness  matrix  before  step  5,  and  follow  the  remaining  steps. 
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III.  RESULTS  AND  DISCUSSION 


In  this  study,  we  considered  three  kinds  of  composites,  Graphite/Epoxy, 
Glass/Epoxy  and  Boron/Epoxy.  The  properties  of  these  composites  are  listed  in  Table  1. 
The  relationship  between  the  transverse  elastic  modulus  {E)  of  each  composite  and  the 
fiber  volume  fraction  {Vf)  for  these  three  composites  is  shown  in  Figure  5.  The  fiber 
transverse  elastic  moduli  of  Glass  and  Boron  are  higher  than  that  of  Graphite,  while  the 
matrix  transverse  elastic  moduli  are  the  opposite.  Therefore,  when  the  volume  fraction  is 
less  than  0.5,  the  transverse  elastic  modulus  of  Graphite/Epoxy  is  higher  than  the  others, 
but  when  the  volume  fraction  is  greater  than  0.5  the  situation  is  reversed.  Figure  6  shows 
the  function  that  best  fit  the  relation  between  transverse  elastic  modulus  and  the  fiber 
volume  flection  in  graphite/epoxy  materials.  The  figures  of  the  other  materials  are  shown 
in  Appendix  C.  (See  Figure  36,Figure  37) 


Table  1.  The  Properties  of  Matrix  and  Fiber  in  the  Graphite/Epoxy,  Glass/Epoxy, 
and  Boron/Epoxy  Composites  Materials. _ 


Property 

Composite 

Graphite/Epoxy 

Giass/Epoxy 

Boron/Epoxy  | 

Epoxy 

Matrix 

Glass  Fiber 

Epoxy 

Matrix 

Boron  Fiber 

Epoxy 

Matrix 

Transverse 

Elastic 

Modulus 

(Gpa) 

15 

5.35 

72.38 

2.75 

400 

2.75 

Poisson’s 

Ratio 

0.49 

0.354 

0.2 

0.35 

0.3 

0.35 

When  the  interface  between  the  matrix  and  fiber  begins  cracking,  the  interface 
will  be  separated  partially,  as  shown  in  Figure  7a.  Figure  7b.  In  order  to  imderstand  the 
stress  distribution  within  the  fiber  when  there  is  a  partial  interface  crack,  a  uniform 
displacement  was  applied  along  the  uncracked  interface  of  the  fiber  as  shown  in  Figure 
7b.  Stress  distributions  for  the  three  fibers  with  50%  or  60%  crack  length  are  shown  in 
Figure  8  to  Figure  13.  The  figures  for  other  crack  lengths  are  shown  in  Appendix  D.  (See 
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Figure  38  through  Figure  43)  From  the  figures,  it  is  found  that  the  stress  distribution  for  a 
half  of  the  crack  length  of  the  fiber  down  to  the  bottom  is  much  lower  when  compared  to 
other  regions.  The  same  results  can  be  shown  in  the  stress  distribution  over  the  both  fiber 
and  matrix  with  50%  of  crack  length  ratio  and  0.36  of  fiber  volume  fraction  for  the  three 
composites  as  seen  in  Figure  14  to  Figure  16.  The  figures  of  other  cases  are  shown  in 
Appendix  E,  Figure  44  through  Figure  49.  Besides  the  part  of  the  fiber  we  mentioned 
above,  the  stress  distribution  at  the  lower  region  of  the  left-hand  side  of  the  matrix  is  also 
negligible.  Therefore,  for  the  analytical  model,  the  lower  stress  regions  were  removed  as 
shown  in  Figure  17,  and  equation  (2-6)  was  used  to  calculate  the  average  composite 
stresses  and  transverse  elastic  modulus  for  a  different  crack  length. 

The  results  of  the  analytical  model  and  the  detailed  FEM  model  with  Vf=  0.36  for 
the  three  composites  are  shown  in  Figure  18  to  Figure  20.  The  results  for  other  fiber 
volume  fractions  are  shown  in  Appendix  F,  Figure  50  through  Figure  55.  They  compare 
the  reduced  transverse  elastic  modulus  ratio  (Ec/Eq)  for  different  crack  length  ratios  and 
different  fiber  volume  fractions.  According  to  the  results,  when  the  crack  length  ratio 
increases,  the  transverse  elastic  modulus  decreases,  and  the  error  between  the  FEM  model 
and  the  analytical  model  for  the  Graphite/Epoxy  composite  is  larger  than  those  for  the 
other  two  composites.  This  is  a  result  of  the  stress  ratio  of  the  region  which  is  neglected 
in  the  analytical  model.  The  stress  ratio  at  these  regions  for  the  Graphite/Epoxy 
composite  is  higher  than  those  for  the  other  two  composites.  So  the  error  between  the 
FEM  and  analytical  models  for  the  Graphite/Epoxy  composite  is  somewhat  larger. 
Nevertheless,  these  two  models  produce  very  compatible  results  and  the  analytical  model 
can  be  used  to  predict  the  effective  transverse  elastic  modulus  for  different  fibrous 
composites  with  a  partial  interface  crack  length. 

As  mentioned  in  Chapter  II,  we  model  springs  in  the  connection  nodes  between 
the  fiber  and  the  matrix  by  adding  the  spring  stiffness  matrix  into  the  system  global 
matrix,  as  shown  in  Figure  21.  For  this  model,  the  spring  constants  kx  and  ky  are  chosen 
from  5  X  10‘°  to  10^.  When  the  spring  constant  is  equal  to  or  larger  than  5  x  lO'®,  it’s 
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assumed  that  there  is  no  interface  crack  between  the  fiber  and  matrix,  and  when  the 
spring  constant  is  less  than  10^,  the  complete  debonding  between  the  fiber  and  matrix 
would  take  place.  Otherwise  there  is  a  partial  debonding. 

Based  on  the  above  concept,  we  define  two  crack  cases,  low  crack  case  and  upper 
crack  case,  as  shown  in  Figure  22.  When  there  is  no  crack  at  point  “A”  and  point  “B” 
starts  to  crack,  ka  is  equal  to  5  x  10’°  and  kb  changes  from  5  x  lO'*^  to  10^.  This  is  the 
low  crack  case.  On  the  other  hand,  when  point  “B”  totally  debonds,  this  means  kb  is  equal 
to  10^,  and  ka  changes  from  5  x  10'°  to  10^.  This  is  the  upper  crack  case.  The  results  of 
these  two  cases  for  different  materials  with  0.36  of  the  fiber  volume  fraction  are  shown  in 
Figure  23  through  Figure  28.  Appendix  G,  Figure  56  through  Figure  67  shows  the  other 
fiber  volume  fractions.  These  figures  illustrate  that  the  spring  constant  and  the  crack 
length  ratio  form  an  exponential  relationship.  The  best  fitting  function  for  each 
relationship  is  also  shown  in  the  figures. 

When  the  spring  constant  decreases,  no  matter  what  the  case  is,  upper  or  lower, 
the  partial  crack  ratio  increases,  and  the  elastic  modulus  will  decrease.  Figure  29  through 
Figure  31  show  the  relationship  of  the  transverse  elastic  modulus  ratio  and  the  crack 
length  ratio  for  the  three  models.  The  other  results  for  different  fiber  volume  fractions  are 
shown  in  Appendix  H,  Figure  68  through  Figure  73.  According  to  these  results,  the 
simplified  FEM  micromechanical  model  with  a  proper  spring  constant  can  also  predict 
the  effective  elastic  modulus  where  there  is  a  partial  interface  crack. 

In  the  above  results,  it  is  assumed  the  shape  of  the  fiber  is  rectangular.  If  we 
change  the  shape  of  the  fiber  to  be  octagon  instead,  then  we  consider  a  new 
micromechanical  model.  Figure  32  shows  the  two  different  micromechanical  models. 
We  calculate  the  effective  transverse  elastic  modulii  of  the  new  micromechanical  model 
with  the  same  material  properties  that  used  in  the  first  simplified  model,  and  compared 
the  results  with  the  first  model. 

Table  2  shows  the  effective  modulii  (undamage)  of  these  two  simplified  models 
for  the  three  different  composites  in  which  fiber  volume  fraction  is  0.36,  and  Figure  33 
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through  Figure  35  shows  the  effective  modulii  of  the  new  micromechanical  model  in 
cracking  condition. 


Table  2.  Effective  Modulis  (Undamage). 


Modulis 

Model 

Rectangular  Fiber  Shape 

Octagon  Fiber  Shape 

Transverse  Elastic  Modulus 
of  Graphite/Epoxy 

7.638  X  10^ 

7.503  X  10^ 

Transverse  Elastic  Modulus 
of  Glass/Epoxy 

6.098  X  10^ 

5.2147  X  10^ 

Transverse  Elastic  Modulus 
of  Boron/Epoxy 

6.538  X  10^ 

5.5310  X  10^ 

Fiber  Volume  Fraction 

0.36 

0.365  1 

20 


2 


a  -  Fiber  Subcell 

b,  c  and  d  -  Matrix  Subcells 


A  Partial  Crack  Situation  of  Subcells 


(b) 


Fiber  Part  of  Subcell  in  Cracking 


Figure  7.  Partially  Separated. 
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Stress 


Figure  8.  Stresses  Distribution  Of  Graphite  Fiber  With  50%  Crack  Ungth. 
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Stress 


Figure  9.  Stresses  Distribution  Of  Graphite  Fiber  With  60%  Crack  Length. 
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Stress 


Figure  10.  Stresses  Distribution  Of  Glass  Fiber  With  50%  Craek  Length. 
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Stress 


Figure  11.  Stresses  Distribution  Of  Glass  Fiber  With  60%  Crack  Length. 
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Stress 


Stress 


Figure  13.  Stresses  Distribution  Of  Boron  Fiber  With  60%  Crack  Length. 
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Stress 


Figure  14.  Stresses  Distribution  Of  Graphite/Epoxy  With  50%  Crack  Length.  (  Vf 

0.36) 
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Stress 


Figure  15.  Stresses  Distribution  Of  Glass/Epoxy  With  50%  Crack  Length.  (F/  - 

0.36) 
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Stress 


xIO® 


Figure  16.  Stresses  Distribution  Of  Boron/Epoxy  With  50%  Crack  Length.  (Vf  — 

0.36) 
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Transverse  elastic  modulus  ratio:  Ec/Eo 


o 


0.8  h 


0.7  h 


0.6 


0.5 


— " :  analy.model 
“o"  :  F.E.M. 


V  o 


-0  0.1  0.2  0.3  0.4  0.5  0.6  0.7  0.8  0.9 

Crack  length  ratio ;  C 


Figure  18.  Transverse  Elastic  Modulis  Ratio  Of  A  Graphiteffipoiy  Composite.  (V> 

0.36) 
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o" :  analy.model 

;  fitting  function  kb=5.6319e10*exp(-9.7324*C) 


Figure  23.  Spring  Constant  Of  A  Graphite/Epoxy  Composite  With  Vf  =  0.36  (Low 

Crack  Case). 
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Spring  const  ka 


Crack  length  ratio  C 


Figure  24.  Spring  Constant  Of  A  Graphite/Epoxy  Composite  With  Vf  =  0.36  (Upper 

Crack  Case). 
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Spring  constant  kb 


Crack  length  ratio  C 


Figure  25,  Spring  Constant  Of  A  Glass/Epoxy  Composite  With  Vj-  =  0.36  (Low 

Crack  Case). 
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Spring  const  ka 


Figure  26.  Spring  Constant  Of  A  Glass/Epoxy  Composite  With  Vj  =  0.36  (Upper 

Crack  Case). 
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Spring  constant  kb 


Figure  27.  Spring  Constant  Of  A  Boron/Epoxy  Composite  With  Vf  -  0.36  (Low 

Crack  Case) 
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Spring  const  ka 


0.7  0.8  0.9 

Crack  length  ratio  C 


Figure  28.  Spring  Constant  Of  A  Boron/Epoxy  Composite  With  Vf  -  0.36  (Upper 

Crack  Case). 


Transverse  elastic  modulus  ratio  Ec/Eo 


Figure  29.Transverse  Elastic  Modulis  Of  Graphic/Epoxy  Composite  {Vf  -0.36). 
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Transverse  elastic  modulus  ratio  Ec/Eo 


Figure  31.  Transverse  Elastic  Modulis  Of  A  Boron/Epoxy  Composite  (V^  -  0.36). 
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Transverse  elastic  modulus  ratio  Ec/Eo 


Transverse  elastic  modulus  ratio  Ec/Eo 


Figure  35.  Transverse  Elastic  Modulis  Of  Boron/Epoxy  For  New  Model  (Vj-  = 

0.365). 
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IV.  CONCLUSIONS  AND  RECOMMENDATIONS 


Both  modified  analytical  and  simplified  FEM  micromechanical  models 
developed  in  this  study  can  predict  the  effective  elastic  modulii  of  composite  materials 
with  a  partial  interface  crack  between  the  fiber  and  matrix.  As  a  result  of  this  study,  the 
following  conclusions  and  recommendations  are  made. 

A.  CONCLUSIONS 


•  A  simple,  analytical,  micromechanical  model  that  neglects  the  non-effective 
load-carrying  zone,  can  predict  the  effective  elastic  modulii  accurately. 

•  A  simplified  finite  element  micromechanical  model  that  contains  interface 
springs  can  also  predict  the  effective  elastic  modulii  accurately. 

•  Different  fiber  shape  micromechanical  models  which  have  the  same  fiber 
volume  fraction,  and  the  same  material  properties,  yield  comparable  results. 

•  When  the  fiber  elastic  modulus  is  much  larger  than  the  matrix  elastic  modulus, 
the  results  of  the  simplified  models  are  more  accurate  than  those  for  the  case 
where  the  fiber  and  matrix  elastic  modulii  are  not  much  different. 

B.  RECOMMENDATIONS 

•  The  micromechanics  models  need  to  be  coupled  with  an  analysis  computer 
program  for  general  composite  structures. 

•  Non-elastic  deformations  of  composites  need  to  be  included  in  the  present 
models. 

•  A  follow-up  research  from  this  study  would  be  to  develop  an  algorithm  for 
determining  spring  constants,  dependend  on  the  state  of  the  interface  crack. 
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APPENDIX  A.  SUBCELL  STRAINS  EQUATIONS 
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♦Equations  10  through  12  <3  =  yjv^ ,  for  a  subscript  “a,b,c,d”  are  subcells 


APPENDIX  B.  MATRIX  [B] 
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APPENDIX  C.  TRANSVERSE  ELASTIC  MODULUS  AND  BEST  FITTING 

FUNCTION. 


X  10 


Transverse  elastic  modulus 


Figure  37.  Transverse  Elastic  Modulus  Boron/Epoxy  and  Best  Fitting  Function 
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APPENDIX  D.  STRESSES  DISTRIBUTION  OF  A  FIBER. 


Figure  38.  Stresses  Distribution  of  a  Graphite  Fiber  With  70%  Crack  Length 


Stress 


Figure  39.  Stresses  Distribution  of  a  Graphite  Fiber  With 
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Stress 


Figure  40.  Stresses  Distribution  of  a  Glass  Fiber  With  70%  Crack  Length 
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Stress 


Figure  42.  Stresses  Distribution  of  a  Boron  Fiber  With  70%  Crack  Length 
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Figure  43.  Stresses  Distribution  of  a  Boron  Fiber  With  80%  Crack  Length 
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APPENDIX  E.  STRESSES  DISTRIBUTION  OF  A  COMPOSITE. 


Figure  44.  Stresses  Distribution  of  a  Graphite/Epoxy  Composite  With  50%  Crack 

Length  (Iy-=  0.49) 


67 


Figure  45.  Stresses  Distribution  of  a  Graphite/Epoxy  Composite  With  50%  Crack 

Length  (k/=  0.64) 
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Figure  46.  Stresses  Distribution  of  a  Glass/Epoxy  Composite  With  50%  Crack 

Length  0.49) 
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Figure  47.  Stresses  Distribution  of  a  Glass/Epoxy  Composite  With  50%  Crack 

Length  (f/=  0.64) 
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Figure  48.  Stresses  Distribution  of  a  Boron/Epoxy  Composite  With  50%  Crack 

Length  (f/=  0.49) 
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Stress 


Figure  49.  Str  esses  Distribution  of  a  Boron/Epoxy  Composite  With  50%  Crack 

Length  (Fy=  0.64) 
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APPENDIX  F.  TRANSVERSE  ELASTIC  MODULUS  RATIO. 


Transverse  elastic  modulus  ratio:  Ec/Eo 


Figure  51.  Transverse  Elastic  Modulus  Ratio  of  a  Graphite/Epoxy  Composite  With 

^^=0.64 
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Transverse  elastic  modulus  ratio:  Ec/Eo 


Crack  iength  ratio  :  C 


Figure  52.  Transverse  Elastic  Modulus  Ratio  of  a  Glass/Epoxy  Composite  With  Vf 

=0.49 
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ulus  ratio:  Ec/Eo 


Crack  length  ratio  :  C 


Figure  53.  Transverse  Elastic  Modulus  Ratio  of  a  Glass/Epoxy  Composite  With 

=0.64 
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modulus  ratio:  Ec/Eo 


ulus  ratio;  Ec/Eo 


Spring  constant  kb 


APPENDIX  G.  SPRING  CONSTANT  OF  A  COMPOSITE. 


Figure  56.  Spring  Constant  of  a  Graphite  /Epoxy  Composite  with  I^=0.49  (Low 

crack  case) 
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Spring  constant  kb 


Figure  58.  Spring  Constant  of  a  Glass  /Epoxy  Composite  with  Fjr-0.49  (Low  crack 

case) 
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Spring  const  ka 


Crack  length  ratio  C 


Figure  62.  Spring  Constant  of  a  Graphite  /Epoxy  Composite  with  i^=0.49  (Upper 

crack  case) 
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Spring  const  ka 


Figure  63.  Spring  Constant  of  a  Graphite  /Epoxy  Composite  with  i^=0.64  (Upper 

crack  case) 
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Spring  const  ka 


Crack  length  ratio  C 


Figure  64.  Spring  Constant  of  a  Glass  /Epoxy  Composite  with  1^=0.49  (Upper  crack 

case) 


87 


Spring  const  ka 


Figure  66.  Spring  Constant  of  a  Boron  /Epoxy  Composite  with  f^=0.49  (Upper 

crack  case) 
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Spring  const  ka 
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Figure  67.  Spring  Constant  of  a  Boron  /Epoxy  Composite  with  J^-0.64  (Upper 

crack  case) 
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Transverse  elastic  modulus  ratio  Ec/Eo 


APPENDIX  H.  TRANSVERSE  ELASTIC  MODULUS  OF  COMPOSITE. 


Crack  length  ratio  C 


Figure  68.  Transverse  Elastic  Modulus  of  a  Graphite/Epoxy  Composite  of  Three 

Models  {Vf  =  0.49) 
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Figure  70.  Transverse  Elastic  Modulus  of  a  Glass/Epoxy  Composite  of  Three  Models 

{Vf  =  0.49) 


Transverse  elastic  modulus  ratio  Ec/Eo 


Figure  71.  Transverse  Elastic  Modulus  of  a  Glass/Epoxy  Composite  of  Three  Models 

(K/  =  0.64) 
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Transverse  elastic  modulus  ratio  Ec/Eo 


Figure  72.  Transverse  Elastic  Modulus  of  a  Boron/Epoxy  Composite  of  Three 

Models  {Vf  =  0.49) 
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Figure  73.  Transverse  Elastic  Modulus  of  a  Boron/Epoxy  Composite  of  Three 

Models  (F.  =0.64) 
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